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I. INTRODUCTION 

Interest in the classical spin systems [ Q] on a two-dimensional lattice has been revived recently by the 

results of a MonteCarlo simulation [ [^] interpreted as evidence of a second order "topological" phase transition, taking 
place for values of the spin dimensionality n > 3. This is unexpected according to renormalization group ideas. Indeed 
the RP"~^ models have the same formal continuum limit as the conventional 0(n) symmetric n— vector spin models, 
therefore they should belong to the same universality class and should not behave differently for n > 3 ( when n = 2 
the i?P"~^ model trivially reduces to the n— vector model). However the global topologies of the spin manifolds: the 
hypersphere S"~^ with antipodal points identified in the case of the RP^~^ model and simply 5'""^ in the case of 
the n-vector model, are different and it has long been known that this might be a reason for different phase diagrams[ 




MonteCarlo studies of these systems, mainly in the n — i case, sometimes with conflicting or not completely 
convincing results are by now numerous [ and they have been augmented by recent more extensive 

simulations on large lattices [ |^,|ll|,|l2|] using cluster algorithms [ jl^] in order to reduce the critical slowing-down. 
On the other hand, high temperature expansion (HTE) studies are still practically absent, the only exceptions being, 
to the best of our knowledge, a series through order 0^ for the internal energy and a series for the mass gap through 
order (3^ in the n = 3 case 

[ §]■ These expansions have been helpful for a first check of MonteCarlo simulation codes, 
and series for other quantities and for other values of n would be equally welcome. 

We have extended to every value of the spin dimensionality n through order /3® the computation of the internal 
energy , and for the first time we have computed series for the susceptibility and the second correlation moment. 

These series are probably not long enough to provide, by their own, convincing evidence about the existence, the 
location and the nature of a possible critical point, but we believe it is useful to make them promptly available so that 
they can serve not only to check MonteCarlo data, but also for future more extensive high temperature calculations. 

We shall explain later why our computational method, based on the Schwinger-Dyson recursion equations[ [Q], 
although very transparent, becomes rapidly cumbersome and therefore is unable, in its present form, to produce 
substantially longer series. 

II. THE HIGH TEMPERATURE SERIES 

Let us briefly describe the model and fix our conventions. 
The partition function of the model is 

Z= fW ds{x)6{s{xf - l)exp[^ E (^(^) • + ^'^))'] (1) 
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The variables of the model are n-component classical spins s{x) of unit length associated to each site x = xici + X2e2 
of a 2-dimensional square lattice, ei and 62 are the two elementary lattice vectors. 

The Hamiltonian and the integration measure have a global 0{n)/Z2 and a local Z2 invariance. Since in two 
dimensions continuous symmetries are unbroken[ [^ ], the most general correlation function < 4'{C) > can be written 
as 

< (j){C) >=< (j){xi,X2,-;Xn;{h,j}) >=< Yl i'^i^^) ■ ^i^j))''" > (2) 

l<i<j<n 

with integer bi^j > 0. The local invariance under Z2, which also cannot break[ p^], implies the further restriction 
that each s{xi) has to appear in (j){C) an even number of times. 

The correlation function (^) may be represented graphically as follows: the lattice points xi, X2, Xn are taken 
as vertices and a line connecting the vertices Xi and Xj is associated to each factor s{xi) ■ s{xj) in (j){C). In terms 
of graphs the local Z2 invariance requires that the degree of each vertex be even. Thus, for instance, the correlation 
< s{xi) ■ s{x2) > vanishes trivially. 

The fundamental two-spin correlation is then G{x2 — xi\ /3, n) =< (s(xi) • s(x2))^ >. 

In particular we have —G{ei;(3, n) = E, the energy per site. 

We also have computed the moments m^'-' {f3, n) of the connected correlations 

C{X2 -xi;p,n)^J2< (^s^{xi)s\xi)- < s''{xi)s'>{xi) > ) (s"(a;2)s''(x2)- < s''ix2)s\x2) > ) G{x2 ~xi;p,n)~ l/r 

a,b 

(3) 



which are defined as follows 



'(/3,n)=^|x|'C(x;/3,n)=^aW/3'- (4) 

X r 

The HTE coefficients for G(ei;/3,n) = Y.r 9r{n)(3'' are: 



go{n) = - 
n 



n-1 

gi[n) = — 



n2(n + 2) 

(n-l)(Ti-2) 
n3(n + 2)(n + 4) 

{n - 1)(72 + I8n - \\r? - + n^) 
n4(n + 2)3(r7, + 4)(n + 6) 

in - l)(n - 2)(528 + 130n - Ylr? - Sn^ 



52 W 

54(71) = 



+ 2)3(n + 4)(n + 6)(n + 8) 



g^{n) = {{n - 1)(284160 + 130496n - 104032n2 - 53344^''' + 688871"* + 5496n^ + 474n^ - 56n^ - 2n' 

/(n^(n + 2f{n + Af{n + 6)(n + 8)(7i + 10)) 

gein) = {n-l){n- 2)(11704320 + 8093952n - 1233088n2 - 1863104^^ - 20 7 7671^ + 103840n^ 
+2621071*^ + 138671^ - lOOn^ - 2n^ + n^")/{n'^{n + 2f{n + ifin + 6){n + 8)(n + 10)(n + 12)) 

grin) = {n - 1)(341118812160 + 428301582336n + 17644511232n2 - 19154965708871^ - 76694446080n'' 
+ 1727682624071^ + 16424658272nS + 192669780871^ - 951227456n^ - 295105184n^ 
-550562671^° + 10001781n" + 1876337ni2 + ISmin^^ + 1527n*^ - 171n^^ + 9n^^ + 7i") 

/{n^{n + 2f{n + Af{n + 6f{n + 8){n + 10)(n + I2){n + 14)) 
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gs{n) = {n~ l)(n - 2)(1271577968640 + 1237547925504n - 87404783616n^ - 44139305932871^ 
-10708273932871^^ + 3754625648071'^ + 17834481104n^ + 440575008n^ - 777645296n** 
-105547274n^ + 1013485371^° + 3591697n" + 31689177^^ + 774977^^ - 28977^'' - 77^^ + 77^^) 
/(n^(77 + 2)^(77 + 4)3(71 + 6)3(n + 8)(7i + 10)(n + 12)(77 + 14)(77 + 16)) 

For 71 = 3 we have (compare with Ref. [ 0]): 

r(e -B 0.1 - ^ , 2 .^2 I 2 34 , 13402 10702 

^iei,p,o; g + ^^p + g^^p + ^g^^p + 467775^ 2280403125^ 47888465625'^ 

12179386 ^7 33996598 ^„ 

-/3 + .o-.o.iino.nir.^o. /3 +■■■ 



142468185234375' 4872411935015625' 



The HTE coefficients for rrS'^\j3,n)^ also called the susceptibility, are: 



n'^(2 + 7i) 



(0) _ 4(77- l)(8 + 377 + 71^) 

"2 773(2 + 77)2(4 + 71) 



(0), , 4(77- l)(96 + 647i + 32772 + 57i3 + 7i4) 

"3 " 7l4(2+7l)3(4 + 7l)(6 + 7l) 



(0)^^^ _ 4(71 - 1)(1 + 77) (3456 + 1968n + 570^^ + 8977^ + 971-* + rt') 
7i4(2 + 7i)4(4 + 7i)2(6 + ?i)(8 + 71) 



ar(7i) 



4^ (71) = (4(n - 1) (122880 + 101888n + 4064071^ + 425287i3 + 356967i^ + 11094n^ 
+ 180771*^ + 16271^ + 1077^ + n^))/(?i''(2 + 77)^(4 + 77)^(6 + ?i)(8 + 77)(10 + ?i)) 



4°^ (77) = (4(71 - 1)(-115015680- 7933132877+ 7460966477^ + 9677286477^ + 440060807t4 

+ 1570220871^ + 751331271*^ + 286201671^ + 64 8 5 6071^ + 8717871^ + 7QASn^° + 3647i" + Wn^^ + n^^)) 

I „^6/^ I \^(a I „^2 



/(77'(2 + 77)''(4 + 7l)-'(6 + 71)^(8 + 77)(10 + 7l)(12 + 77)) 



4°^ (77) = (4(71- l)(43104337920 + 438669803527t- 540706406471^- 1500234547277^ + 376586752071^ 
+887809792077^ + 428230528071^ + 11968429127t^ + 326380672n^ + 97376320n^ 
+2212316871^° + 322842277" + 29247271^2 + 160587i" + 56677^^ + 2377^^ + n^^)) 

7(77^(2 + 7l)^(4 + 7l)3(6 + 77)3(8 + 7l)(10 + n)(12 + 77)(14 + 7l)) 
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4°)(n) = (4(n - 1)(-94746307461120 - 126660674322432n - 5226623926272n2 + 66792515567616n^ 
+32795171340288n'' + 481 94863 104n''^ - 862398014464n^ + 2921457334912n^ 
+2239266005664n^ + 7907584401 12n^ + 18502955169671^° + 37818452512n" + 777687597071^^ 
+139524797171^3 + 1845880287i" + 16678488n^^ + 985722n^^ + 3665071^^ + 952n^^ + 32n^^ + 71^°)) 

7(71^(2 + 7i)^(4 + 7i)*(6 + nf{8 + nf{10 + n){12 + n){U + n)(16 + n)) 

For n = 3 these formulae give: 

m i/^,3j - g + 45/^+4725'^ + 70875'^^ + 5457375^^ + 15962821875^ 
37545856 . 10273872032 7 934133719808 g 
''"335219259375^ 427404555703125^ 183909666174609375^ 

The HTE coefficients for m^'^\f3,n), the second correlation moment, are: 

7l2(2 + 7l) 



(2), . _ 4(71- l)(28 + 87l + 7l^) 
„3(2 + n)2(4 + 7l) 



(2), , 4(71 -l)(624 + 344n+ 1247i2 + 1571^ + 71^) 

«3 W = 



71^(2 + 7l)3(4 + 7l)(6 + 7l) 

A(n - 1 Vfi999zL 4- .t^7«.t^fir7 -I- ?t77finTi2 _|_ 1 : 



(2) , , _ 4(71 - 1)(52224 + 5785671 + 3776071^ + 1320071^ + 284471* + 37671^ + 2571^ + ti^) 



n5(2 + 7i)4(4 + 7i)2(6 + 7i)(8 + n) 

a^^^ (n) = (4(71 - 1) (1044480 + 155340871 + 139417671^ + 69276871^ + 237584n* + 50998n^ 
+710771*^ + 64271^ + 3071^ + 7i^))/(7i*^(2 + 7i)^(4 + 7i)^(6 + 7i)(8 + 7i)(10 + n)) 



a^^\n) = (4(71 - 1)(420249600+ 109246464071+ 152355532871^ + 125585561671^ + 7104977287z'^ 
+28557206471^ + 8426452871*^ + 1840414471^ + 290909271^ + 32705471^ + 2580371^° + 13547i" + 4471^^ + 71^^)) 

7(71^(2 + 7i)^(4 + 7i)3(6 + 7i)^(8 + 7i)(10 + 7i)(12 + n)) 



a^^\n) = (4(71- 1) (83979141120 +20430943027271 +29172820377671^ + 30010984857671^ 
+23251085465671* + 13105306240071^ + 5480046937671*' + 1732433724871^ + 421261801671^ 
+79477208071^ + 11486543271^° + 12396858ri" + 97 7 5 3272^^ + 5502871^3 _^ 2 1067^" + 5371^^ + n^^)) 

7(71^(2 + 7i)^(4 + 7i)3(6 + 7i)3(8 + 7i)(10 + 7i)(12 + 7i)(14 + n)) 



ai\n) = (4(71- 1)(-58788371496960+ 836470466150471+ 28017111859200071^ + 47950452051148871^ 

+46420291141632071* + 329227765829632tj;"^ + 188047485044736r?.^ + 8702963842406471^ 
+3217373044352071^ + 945655868582471^ + 22198000183687ii" + 4193256525767i" + 6393045419271^2 
+7800108 77671^ + 74699521271^* + 5484162071*'^ + 301399271**^ + 12045471*^ + 33747ii* + 677ii° + ti^O)) 

7(71^(2 + 7i)^(4 + 7i)*(6 + nf{S + 7i)2(10 + 7i)(12 + 7z)(14 + 7i)(16 + n)) 
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In particular, the HT expansion of m'^'^\j3, 3) is: 



(2\,n " ^ -ioo n2 2896 1123712 ^4 67018144 2023066384 

3^ = B A B A B A B A B A B 

VH, J 45^^^ 4725'^ 70875^ 81860625^ 15962821875'^ 1676096296875^ 

12824336768 ^7 18110407484144 ^„ 



38854959609375' 208430954997890625' 

A correlation length may be defined, as usual, in terms of the ratio of m*^^'(/3, n) and m(°^(/3, n). 

Let us notice that a few simple checks of the formulae are possible: all HTE of the connected correlations have 
to vanish for n = 1 because of the triviality of the model. For n — 2, the expansions should reduce to the 

corresponding ones for the 0(2) (or XY-) vector model. Finally, for rt = 3 the HTE of C(ei;/3, n) agrees with the 
calculation of Ref. [ 0]. 

Our HTE have been computed from the Schwinger-Dyson equations of the model, an infinite system of linear 
equations among the correlation functions. The generic equation, which may be deduced following closely Ref. [ JT^], 
has the structure 

n 

< '^(^) >= -9 te(< ^(^P >-< ^(^P >) + (^12 - 1) < </'(Cl2,12) > - < ^(Ci^ly) > 

/ 1 n~ til ^ L . _ 



" + 51 - „ 



(5) 



Here we have assumed that the vertices xi and X2 are connected by one line at least, gi is the degree of the vertex 
Xi , bij the number of lines connecting the vertices Xi and Xj , < ) > denotes the correlation function obtained 

from < 4>{C) > by removing a factor 5(2:1) • s{x2) and replacing it by s{xi) ■ s(a;i+p)s(x2) • s(a;i+^), namely 

Mn-\ j,/r<\ ^(^'^'> • 5(^:1+^)^(^2) ■ s{xi+^) 
'^(^-^^'^^^^ six,) ■ six2) 

and analogously 

</.(C+)=</.(C)((s(xi).s(xi+^))2 
(s(xi) • s(a;2))2 



The HTE of the correlation < (j){C) > is obtained solving iteratively eqs.(y) by the same procedure as in the case of 
the 71— vector model[ [|4j]. Here however, a difficulty is met: while in the case of the n— vector model a large fraction 
of the graphs generated after the first few iterations can be neglected, in this case, due to the local Z2 symmetry, 
all graphs contribute nontrivially to the final results and therefore must be recorded. Thus the required computer 
memory rapidly becomes exceedingly large and it is difficult to push the expansion beyond the 8-th order. However 
not all the blame should be laid upon the computational technique since the combinatorial complexity of the expansion 
is really higher and of a faster growth with the order than in the n— vector case. It is also interesting to recall that 
analogous difficulties were met when performing strong coupling expansions in the Hamiltonian formalism[ |^]. 

A simple analysis of the series by ratio and Fade approximants methods [ [|l^] (see Fig.l and Fig. 2 ) suggests 
the existence of a critical point when n « 2, but, unfortunately, the series seem to be not long enough to warrant 
any reasonably safe conclusion when n = 3 or greater. To be sure, for various values of n there are some Fade 
approximants of the susceptibility having a real positive singularity or a complex conjugate pair of singularities 
nearby the real positive fB axis and in the expected position. The same happens for the logarithmic derivative of the 
susceptibility. These poles however, at this order of approximation, are not stable enough to enable us to exclude 
the possibility of an artifact of low order approximants to mimic the steep increase of the susceptibility. Thus some 
completely different scenarios are still compatible with our series, for instance: 

a) in analogy with the behavior of the n— vector model[ |l^] a critical point exist for n ~ 2. As n is increased and 
varied through some n < 3, the critical point might split into an unphysical pair of complex conjugate singularities 
so that the model becomes asymptotically free for n > 3. This conjecture might be supported both by the alternate 
ratios plots of Fig. 1, which seem to show the onset of an oscillatory trend [ ||l^] and by some Fade approximants to 
the susceptibility or its log-derivative whose nearest singularities in the right half fB plane are complex. 

b) a critical point exists for all n as sugg ested by Ref.[ §]. 



5 



ACKNOWLEDGMENTS 



Our thanks are due to Sergio Caracciolo for suggesting to undertake this computation and to Alan Sokal for 
further encouragement and useful discussions. We also are indebted to U. Wolff for a useful discussion and for kindly 
permitting us to use his unpublished MonteCarlo data in Fig. 2. Finally we thank A. J. Guttmann and G. Marchesini 
for carefully reading a draft of this note. Our work has been partially supported by MURST. 



[1] W. Maier, A. Saupe, Z. Naturforsch. 14A, 882 (1959); 

15A, 287 (1960); P. A. Lebwohl, G. Lasher, Phys. Rev. 

A 6, 426 (1972); M. J. Stephen, J. P. Straley, Rev. Mod. 

Phys. 46, 617 (1974); P. G. de Gennes, The Physics of 

liquid crystals , (London, Oxford University Press 1974). 

M. Plischke, M. Bergersen, Equilibrium Statistical Physics, 

( Englewood Cliffs, Prentice Hall 1989). 
[2] H.Kunz and G. Zumbach, J. Phys. A 22, L1043 (1989); 

H.Kunz and G. Zumbach, Phys. Lett. B 257, 299 (1991). 
[3] J. Kogut, M. Snow, M. Stone, Nucl. Phys. B 215[FS7], 

45 (1983). 

[4] R. D. Mountain, Th. W. Ruijgrok, Physica 89A, 522 
(1977). 

[5] S. Duanc, M. B. Green, Phys. Lett. B 103, 359 (1981). 

[6] S. Solomon, Phys. Lett. B 100, 492 (1981). 

[7] M. Fukugita, M. Kobayashi, M. Okawa, Y.Oyanagi, A. 

Ukawa, Phys. Lett. B 109, 209 (1982). 
[8] D. K. Sinclah, Nucl. Phys. B 205[FS5], 173 (1982). 
[9] M. Caselle, F. Gliozzi, Phys. Lett. B 147, 132 (1984); M. 

Caselle, Zeit. fur Phys. C 28, 233 (1985). 
[10] C. Chiccoli, P. Pasini, C. Zannoni, Physica 148A, 298 

(1988). 

[11] U. Wolff, unpublished work 

[12] S. Caracciolo, A. Pelissetto, A. D. Sokal, work in progress. 
[13] J. Goodman, A. D. Sokal, Phys. Rev. Lett. 56, 1015 (1986); 

R. H. Swendsen, J. S. Wang, Phys. Rev. Lett. 58, 86 

(1987); U. Wolff, Phys. Rev. Lett. 62, 361 (1989); J. 

Goodman, A. D. Sokal, Phys. Rev. D 40, 2035 (1989); A. 

D. Sokal, Nucl. Phys. B (Proc. Suppl.) 20, 55 (1991). 
[14] P. Butera, M. Comi and G. Marchesini, Nucl. Phys. B300 

[FS22], 1 (1988); P. Butera, M. Comi and G. Marchesini, 

Phys. Rev. B 41, 11494 (1990). 
[15] N. M. Mermin, H. Wagner, Phys. Rev. Lett. 17, 1137 

(1966). 

[16] S. Elitzur, Phys. Rev. D 12, 3978 (1975). 

[17] G. A. Baker and D. L. Hunter, Phys. Rev. B7, 3346, 
3377(1973); D. S. Gaunt, A.J. Guttmann in Phase Tran- 
sitions and Critical Phenomena, C. Domb and M. S. 
Green Eds. Vol. 3, (London, Academic Press 1974); A. J. 
Guttmann in Phase Transitions and Critical Phenomena, 
C. Domb and J. Lebowitz Eds. Vol.13, ( London, Academic 
Press 1989). 



FIG. 1. Alternate ratios f,(n) = (a'-"\{n)/a[°lj^{n)y^^ of 
the expansion cocfhcionts of the susceptibility for various val- 
ues of the spin dimensionality n are plotted versus 1/s. Going 
from the lower plot to the upper, we have n = 2, 3, .., 8. 
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FIG. 2. The susceptibility torn = 3 at order vs. /3. The 
continuous curve shows the [4/4] Pade approximant (which 
is singular at /? « 5.175 ± 0.315i). The dashed curve shows 
the sum of the susceptibility series truncated at order /3®. 
The squares represent data from the MonteCarlo simulation 
of Ref. [ [^]. The triangles represent unpublished data from 
a MonteCarlo cluster simulation performed by U. Wolff [ | pT[ ] 
on lattices of size up to 256^. 



